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Chiral polymers are ubiquitous in nature and in the cellular context they are often 
found in association with membranes. Here we show that surface bound polymers 
with an intrinsic twist and anisotropic bending stiffness can exhibit a sharp continu- 
ous phase transition between states with very different effective persistence lengths 
as the binding affinity is increased. Above a critical value for the binding strength, 
determined solely by the torsional modulus and intrinsic twist rate, the filament 
can exist in a zero twist, surface bound state with a homogeneous stiffness. Below 
the critical binding strength, twist walls proliferate and function as weak or floppy 
joints that sharply reduce the effective persistence length that is measurable on long 
lengthscales. The existence of such dramatically different conformational states has 
implications for both biopolymer function in vivo and for experimental observations 
of such filaments in vitro. 

I. INTRODUCTION 

Both prokaryotic and eukaryotic cells possess a rich variety of filamentous proteins with 
a wide range of compositions, sizes and flexibility [I]. These biopolymers play a vital role 
in a number of critical cellular functions including maintaining structural integrity, serving 
as a template for cell growth, cell adhesion and motility, mechanical signal transduction 
and cell division. In many of these roles, the association of the polymers with membranes 
is critical and necessary for function. For example, the actin rich cell cortex in eukaryotic 
cells is a major determinant of the cell's structural integrity and mediates adhesion and 
signal transduction. A number of different actin binding proteins link the cortex with the 
membrane both directly, and indirectly through other proteins and complexes [2j. In the 
axons and dendrites of nerve cells, members of the MAPI family link microtubules, which 
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are critical for vesicle transport and mechanotransduction, to the plasma membrane [3]. In 
red blood cells, the remarkably flexible spectrin-actin network, responsible for cell shape 
regulation, is tightly associated with the membrane by a number of proteins such as ankyrin 
and dematin [4j. Similarly, the structural support lamin networks inside the nucleus are 
attached to the nuclear envelope by LAP family proteins [5]. During cell division, forces 
required for pinching in the division plane are generated by shrinking actin contractile rings 
attached to the membrane |6]. A similar phenomena occurs in bacterial cell division where 
the bacterial filamentous protein FtsZ makes up the contractile ring [7] which can be linked 
to the membrane by FtsA, ZipA and ZapA [SI Another bacterial protein MreB which 
provides structural support and guides cell wall synthesis is attached to the membrane via 
RodZ, MreC/D [9l[Tn] or even directly jUj. While these and several more examples illustrate 
the common nature of filamentous proteins binding to membranes in vivo, it is also known 
that many filamentous proteins can show direct interactions with membrane components in 
vitro typically through electrostatic or hydrophobic interactions jl] . This is of fundamental 
importance when considers that many experimental setups that study such filaments in vitro 
have them in close proximity or immobilized/fixed on a surface. Clearly, understanding 
the association of filamentous proteins with membranes and surfaces is important from a 
scientific and technological perspective. 

To date, there has a been a large amount of work studying polymers at interfaces and sur- 
faces [121 US] • However, many of these studies tend to neglect key features of real biopolymers 
that could have qualitative effects on their structure and their interactions with surfaces. 
Firstly, a vast majority of biopolymers are chiral with an intrinsically helical structure. This 
helicity can manifest itself in a number of fascinating ways such as the cholesteric liquid 
crystal and the blue phases of concentrated DNA solutions and the self-limiting bundle size 
of aggregates of generic chiral biopolymers such as actin [HHUj. Furthermore, classical 
models of polymers such as the worm-like chain model do not take anisotropies in elastic 
moduli or finite widths of the biopolymers into account. A number of recent efforts have 
therefore focused on the statistical mechanics of elastic "ribbons" with a finite aspect ratio 
as a more realistic model for typical biopolymers or other protein aggregates [T81428] . The 
ribbon-like nature of the polymers can lead to interesting effects including layering transi- 
tions in highly anisotropic condensates and the existence of an underlying helical structure 
even in the absence of twist [221 ED]- Thus the introduction of helicity and anisotropy can 
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result in qualitative changes in the structure and dynamics of both single filaments and their 
aggregates. 

Given the ubiquity of membrane and surface adsorbed biopolymers, in this paper, we take 
the first step toward understanding how hclicity and anisotropy can infiucncc the structure 
and dynamics of polymers that are bound to surfaces. Rather than study a polymer with 
more detailed resolution, for example a full atomistic representation, we focus on the generic 
effects of intrinsic helicity on polymer adsorption using a minimal, coarse-grained description 
of molecular structure. We show that the interplay between helicity and surface interactions 
can have profound effects on the structural and thermodynamic properties of the polymer. 
Wc thus probe the consequences of helical structure on the large-length scale properties of 
absorbed chains, which is of potential relevance to a broad class of biomacromolecules. We 
first build a general model that allows us to prescribe the thermodynamics of a helical and 
anisotropic filament binding to a surface of arbitrary curvature. For the rest of the paper, 
we then treat only planar surfaces. We first consider the infinitely stiff limit where the 
anisotropy becomes irrelevant and only the helicity plays a role. 

The key conclusion from this analysis, is that the structure and conformational fluctua- 
tions of a helical polymer arc critically sensitive to the attractive interactions with a surface. 
We find that the filament undergoes a rapid and continuous phase transition from a strongly 
bound untwisted state to a weakly bound state for a critical value of the binding strength.The 
transition occurs via a proliferation of "twist- walls" , regions where the chain unbinds locally 
from the surface and rapidly twists a single turn. We then consider how this transition alters 
the conformations and fiuctuations of an anisotropic fiexible filament on a planar surface. 
The coupling between the twist and bend degrees of freedom introduced by the anisotropy 
and binding results in surface bound conformations with highly heterogeneous bending stiff- 
ness below the critical binding strength. The transition between strongly and weakly bound 
states is then accompanied by a sharp drop in the effective bending rigidity of the polymer. 
We demonstrate these conclusions, first, through an exact, "zero-temperature" analysis of 
the continuum mechanical model of the filament that ignores thermal fiuctuations of twist. 
We then extend our analysis to finite-temperature by way of Monte Carlo simulations, which 
demonstrate that the sharp transition between the strongly-bound and weakly-bound states 
of helical filaments persists in the presence of both bending and torsion fiuctuations of the 
chain backbone. Thus, we predict that strong couphng of the twist degrees of freedom of a 



4 



helical polymer to surface interactions, not only restructures the torsional state of the chain, 
but it also leads to a dramatic transition in the conformational fluctuations of the entire 
chain backbone. This suggests that experimentally measured properties of biopolymers ab- 
sorbed to a surface such as the persistence length and end-to-end size are crucially sensitive 
to the strength of surface interactions. 



This paper is organized as follows. In Section, II A, we build the general model for a heli- 



cal and anisotropic filament binding to a surface of arbitrary curvature and in Section. II B 



we describe our Monte Carlo simulation procedure. In Section II C we analyze the model 
for planar surfaces in the infinitely stiff filament limit. We specifically study the transi- 



tion between the strongly bound and weakly bound states of the filament. In Section IIP 



we study how the transition influences the in-plane bending fluctuations of an anisotropic 
surface adsorbed filament. We conclude with a short discussion in Section HTTl 



II. CONTINUUM THEORY OF SURFACE-BOUND HELICAL FILAMENT 



A. Model 



Here we introduce to a continuum model of a helical filament bound to an attractive 
surface. For the purposes of generality, we consider the filament bound to a cylindrical 
surface of radius, r. This is a natural geometry to describe the binding the cylindrical cell 
wall of a bacterium, as in j9j , though we will specialize in the next section to the case of 
planar surfaces, r — )• oo. 

The filament is modelled by a helical strip, whose cross-section is rectangular with a 
thickness, t and width, w. In the absence of external forces or interactions with the surface, 
the wide axis of the filament has a natural rotation rate, Uq, along the long axis of the filament 
(see Fig. [2]). Both thermal energy and surface interactions distort the filament geometry, 
leading to elastic costs which we described with the Kirchoff-Love of elastic beams [31] . 
The deformations of the filament are described by rotations of the material frame along the 
arc-length, s, of the filament backbone. We choose ei (62) to describe the local orientation 
of the wide (thin) axis of the filament, and 63 = t to describe the tangent. 

Here, we compute the mechanical energy of the filament in terms of an angle, 0{s), 
describing in-plane orientation of t and an angle, ip{s), describing the orientation of the 
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FIG. 1: Diagram defining the coordinates used to describe an intrinsically helical filament, (a) 
A section of an intrinsically twisted filament which is forming a helix in a cylindrical coordinate 
system, where t, ei and e2 are the material frame basis vectors. The grey curve is the center line 
of the filament, (b) The projection of the tangent vector, t, along z and (p defines of the bending 
angle 0{s). (c) The rotation of the material coordinates ei and e2 defines of the twist angle V'('S)- 

wide axis of the filament with respect to the surface normal. If the center line of the 
filament, r(s), is confined to a helical surface (with long axis centered on the x = y = Q 
axis), then the tangent vector can be described in terms of the angle 9{s), where 

S^r = t = cos 6'(s)i + sin 6'(s)0, (1) 

where = z x r(s)/r is the local, azimuthal direction on the cylinder. The curvature of 
the filament center is determined by S^t, which is non-zero when 0' 7^ and the filament 
curves in the surfaces of the cylinder, or when sin^ 7^ 0, and the filament, following a 
straight, or geodesdic, path on the cylinder is deflected normal to the cylinder. That is, the 
direction changes along the arc-length as the curve is convected around the cylinder axis 
as dscj) = — 5i^f where the radial direction is given by f = [r — z{r ■ z)]/r . Hence, 

a,t = -^^f + ^'(t xf), (2) 
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FIG. 2: Schematic of a helical filament with an anisotropic cross section. On the right, we show a 
series of helical configurations for increasing strength of surface interactions, which show progressive 
unwinding of the helical twist. 

where t x f = sin 6z — cos Oip. To describe the local anisotropy of the filament cross-section 
we define the material directions, perpendicular to t, 



ei = sin-?/'f — cos'?/'(t x f) 



(3) 



and 



62 = COS ipf + sm'ijj(i X f). (4) 
Defining Ki = ■ dst as the rate of back bone bending into the direction e^, we have, 

Kl = HSm{lp — S), K2 = Kcos{ip — 6), (5) 

where 



3/\2 



+ 



sin^^ 



(6) 



is the square curvature of the central axis and tan 5 = r9'/ sin^ ^. Finally, we extract the 
rate of cross section twist, u, from the rate of rotation of ei and 62 around the tangent. 



w = 62 ■ dsei = i) , 

Zr 



(7) 
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From these quantities we define the following elastic energy for the helical filament, 



^Jds + + K{u - Wo)'] , (8) 



-^mech 2 

where Ci and C2 are bend moduli for bending around the respective wide and thin axes of 
the filament cross section and K is the twist modulus. Modeling the filament interior as an 
isotropic, elastic medium these three constants are related by 

Mwtf EwH Ewt^ 

^ = 2t 2 I 42^ ' = "To"' ^"^ = "To" 
TT-^lW^ _|_ 12 12 



where E and /i are the Young's and shear moduli of the of the elastic medium |10] . 

Next we introduce an effective potential to describe the thermodynamics of interactions 
between the helical filament and the attractive wall. We assume that filament maintains 
a stronger state of binding when the surface is in close contact with the wide axis of the 
filament, that is, the minimal energy configuration aligns the 62 with the surface normal, r. 
Thus, we consider a simple energy gain per unit length proportional to (62 ■ r)^, 

Sint = j dsV sin^ij), (10) 

where V describes the binding affinity difference (per unit length) between "face on" and 
"edge on" binding of the filament. We note that similar models have been applied to 
study other aspects of biopolymer behavior, such a base-pair stacking interactions between 
bound pairs of DNA [321 |33]. However, to our knowledge, our study is the first to use 
an intrinsically-twisted elastic filament model to study the frustrations between surface- 
interactions and mechanics of helical polymers, and the first to study structural transitions 
driven by these interactions. 

In the remainder of this article, we focus on the case of a filament bound to a planar 
surface (r — )• 00). The lack of surface curvature greatly simplifies the coupling between 9 
and ip, 

E{r ^ 00) = jds (^y cos2 7/^+ y sin^T/^jl^f + y(7/;'-Wo)' + V^sin2v]. (11) 

However, due to the anisotropy of bending response, an important coupling remains be- 
tween in-plane backbone fluctuations and rotations of cross section. It is easier to bend 
the filament about the thin axis, and hence, rotation of ip along the backbone correlates to 
variations stiffness to bending in the plane. We find below this coupling the exists signifi- 
cant consequences for the sensitivity of contour fluctuations of the filament backbone to the 
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state of filament binding. We now describe our Monte Carlo simulations that serve as an 
independent test of our analytic results. 



B. Numerical Simulations 



In this section we describe our numerical setup to study an adsorbed filament on a flat 
surface. In the subsequent section, we will study the analytical solution to the equation of 
mechanical equilibrium for absorbed helical filaments. While this provides an exact descrip- 
tion of the minimal-energy states of the model, it neglects thermally-induced fluctuations of 
the torsion as well as the backbone bending. It is reasonable to expect that at finite tem- 
perature, the presence of both types of fluctuations modify the "zero-temperature" analysis 
of the theory. Thus, we carry out our numerical simulations of the model to quantify the 
contributions of finite-temperature fluctuations absent from the analytic solutions of the 
theory. Monte Carlo simulations were performed to find equilibrium configurations of a 2d 
chain of length L that is parsed into discrete segments or plaquettes. Each plaquette is of 
unit length and in general can undergo bending and torsional deformations. Plaquette i has 
associated with it a pair of angles (6*4, ipi), which represent the local orientation and twist 



of the chain respectively. A discrete version of the energy given by eq. (11) is implemented 
to evaluate the energy of each chain configuration. 



E = -Y^{Ci cos'(V5i,»+i) + C2 sin2(V'i,i+i)) [Mi^i+if + iC(A^i,i+i - u^f + Vsm^iiPi). 

i 

Here i refers to the segment index, refers to the average twist between nearest neigh- 

bor plaquettes, A refers to finite differences between neighbors along the chain and elas- 
tic parameters are in units of fc^T. Configurations of the chain were accepted (rejected) 
upon comparing the trial configuration energy to the previous configuration energy using 
a METROPOLIS algorithm with the usual Boltzmann weight (exp(— AE')) [? ]. Specifi- 
cally we look for configurations z'^^^^\ where t denotes the Monte Carlo step, that represent 
equilibrium configurations of the chain. Trial configurations of the chain are accepted/reject 
by the following criterion. 
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I z' with probability r(2;, z') = min(l, 
\ z if To > r{z^ z') 

Here z is the previous configuration at step t and p{z') is the Boltzmann weight of the trial 
configuration at step t + 1. Also, tq is a random number chosen from a uniform distribution 
between and 1. Updates of either degree of freedom {Oi or ipi) per plaquette were done 
independently with equal probability for each Monte Carlo step. Update acceptance for trail 
configurations of the chain were between 30% and 40% of the total amount of Monte Carlo 
steps (10^ steps per chain). Data was averaged over ~ 10^ equilibrium configurations of the 
chain. Equilibration of the chain was determined by measuring the rotational diffusion of 
the end to end vector using the condition, 

R{t) ■ R{0) = 0, 

where t is the Monte Carlo step. Usual equilibration times were of the order 10^ steps. 
Chain lengths of = 32 were used in all numerical results unless otherwise stated. 



C. Thermodynamics of surface-induced untwisting 

In this section, we analyze the thermodynamics of surface binding of the helical filament, 
initially ignoring the coupling between in-plane orientation and filament twist described in 



eq. (11). This can be justified in the limit of stiff filaments, where stiffness is characterized 



by the 2D persistence length, 

^p-2 , Kli) 

which is the effective persistence length of a planar helical filament twisting at a constant 
rate [41j In the limit L ^ £p, the backbone remains effectively rod-like, and we may assume 
& ~ 0. 

With this approximation in mind, the remaining degrees of freedom describe filament 
twist, and the contact between the wide edge of the filament and surface, the final two 



terms in eq. (11). The underlying frustration of helical filament binding is straightforward to 
understand from these terms. While the twist elastic energy is minimized by a constant rate 
of rotation, = cjqs, the binding energy is minimized by locking-in to constant orientations, 
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if) = nn. The equation of motion describing the minimal energy configurations of is 
well-known to the analyses of ID models of incommensurate solids [M] and rigid pendula, 



l^f = A-2(sin2V^ + e), (13) 
where e is a positive-valued "constant of motion" and 

is a length scale defined by the ratio of twist modulus to strength of surface binding. The 
limit, e — 7- 0, describes "lock in" to an untwisted state of face-on binding, ip = nn. For small, 
but finite, values of e, the solution becomes inhomogeneous. The filament locks into near 
perfect surface registry, with ~ nn, over spans of arc length much greater than A. These 
nearly "commensurate" domains are separated by rapid jumps, twist walls, where ip rapidly 
jumps by tt over an arc-distance of order A. For large e, the filament twists homogeneously. 

The transition from locked-in, untwisted filaments, to helical filaments has a highly non- 
linear, and critical, sensitivity to V and the elastic cost (per unit length) to untwist the 
filament, proportional to Kuq. To analyze this, we follow the analysis of the Emery and 
Bak [35j to derive an effective theory for filament binding in terms of the length-averaged 
rate of twist, {ip') = 2tt /P, where P is the pitch of rotation. This pitch is twice the distance 



between twist walls and can be computed from eq. (13) as function of e 



P{^) = A t!,, = 4Ae-V^ir(-e-^), (15) 

Jo V W) + e 

where K{k) is the complete elliptic function of the first kind. The asymptotic behavior 
near the lock-in transition has the form P(e — 0) ^ — 2Aln(e/16), from which we express 
the constant, as the approximate of function of pitch, e(P) ^ 16e~^/^^. Substituting the 



solution, eq. (13), into eq. (11) and integrating over one helical turn, we may derive an 



expression for the effective elastic energy per unit length as a function of P, 

E{P)/P ~ {SV\ - 2tiKujq)/P + WXe-^/'^^/P + Cq. (16) 

Note that the density of twist walls is 2/P, so that the first term in parentheses is propor- 
tional the energy per twist wall. The cost of the wall corresponds an increase in potential 
energy over a length of order A, while a twist wall relaxes the torsional strain on the un- 
twisted filament, leading the negative contribution to twist wall energy. The second term 



11 



can be thought of as the cost of exponentially screened repulsive interactions between twist 
walls. 

It is straightforward to show that when the elastic energy gain is smaller then the potential 
energy cost of a twist wall, the untwisted, locked-in state (P — )■ oo) is stable. This occurs 
when V > Vc, where 

v;^^. (17) 

However, when V < Vc, the energy per twist wall is negative and the filament begins unbind 
and wind through the incorporation of twist walls. For small — \^|, at the unbinding 
threshold, the equilibrium pitch grows rapidly as 

P^\n\Vc-V\. (18) 

This indicates that the filament undergoes a rapid and continuous phase transition from 
the untwisted and strongly-bound state to the helical and unbound state at V = Vc- The 
dependence of the mean-twist, 2tt/P, can be calculated using results summarized in the 
Appendix. The profile of {ip') = 27r/P vs. V as calculated in absence of thermal fluctuations 
is shown in Fig. |3] (a). In Fig. |2]we show examples of torsional state of a helical filament 
with a rectangular cross section as the unwinding transitions is approached from below. 
Simulation data is plotted on top of Fig. [3] (a) which is in excellent agreement with this 
unwinding transition. 



D. Conformational collapse of desorbing helical filament 

In this section, we consider the effect of the transition between the strongly-bound, un- 
twisted state to the helical state on the backbone fluctuations of the surface confined fil- 
ament. Clearly, when the filament is strongly-bound in the face-on configuration in-plane 
bending, around the wide edge of the filament cross section, is suppressed. As described 
above, just below the binding threshold, a low density of twist walls populate the length 
of the filament. As the wide axis is normal to the surface in these regions, they act as 
localized weak spots for in-plane bending, whose bending stiffness is reduce by a factor of 
{t/wY relative to the face-on spans. Hence, in the weakly-bound, helical state, the bending 
stiffness becomes highly heterogeneous. 
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FIG. 3: Both figures together highhght mechanically what is taking place for both degrees of 
freedom as the filament becomes strongly bound to the substrate. In (a), the solid curve shows 
the predicted dependence of the net filament twist on the strength of surface interaction. Solid 
points with error bars display simulation results at finite temperature. (Red) wjt = 2.0, (Purple) 
w/t = 4.0 and (Blue) w/t = 8.0. In (b), a plot of the angular diffusion, for w/t = 4.0, associated 
in-plane orientational fluctuations (eq. ( |20[ ) of the filament tangent along its backbone for fixed 
profiles of equilibrium twist il^{s) at various levels of surface binding. Here, V'(O) = for each 
configurations. 

This for sufficiently anisotropic filament cross sections, this heterogeneity will have a 
profound impact on the bending fluctuations of the filament backbone in the plane. These 
fiuctuations are observable in the tangent-tangent correlations, 



(t(si) ■ i{s2))^ = exp [ - (|e(si, S2)H^/2] . 



(19) 



where the notation (■)^ denotes an average for a fixed twist profile, il){s). From the mechan- 
ical energy for bending we may derive the following simple result for the "diffusion" of 9 
along the contour length, 

1 



{\e{s,,S2)V)^=2 ds' 



S2 



(20) 



£iCos2V^(s') +^2sin2V^(s')' 
where ii = 2Ci/kBT and £2 = 2C2/kBT, are 2D persistence lengths corresponding to bend- 
ing around, respective, thin and wide directions of the cross section. The experimentally 
measurable tangent-tangent correlation function, C^(s2 — si) = (cos[6'(s2) — 9{si)]), is ob- 
tained by averaging (t(si) ■t(s2))^ over all torsional states. According to the analysis above, 
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FIG. 4: In (a), a plot of the (logarithm) of tangent-tangent correlation function of helical filament 
with bending anisotropy C1/C2 = 16 and wq^i = 1- Three different surface binding strengths 
are shown: (blue) vanishing binding strength; (red) near-threshold binding strength; and (yellow) 
above-threshold binding strength. In (b), a plot of the (logarithm) of the ensemble average tangent- 
tangent correlation along the chain contour n for 3 different binding strengths, here C1/C2 = 64. 
As we can see the oscillations about and average effective persistence length is present in good 
agreement with analytical results. 

distinct configurations of ip{s) at a given binding energy are periodic functions of s (mod- 
ulo shifts by vr), differing only by a uniform translation along the backbone (modulo P/2). 
We may perform the averaging over torsional configurations by fixing the twist angle at a 
reference point, say il){s = 0) = 0, and averaging (t(so) ■ t(so + s))^ over a span P/2 along 
the backbone for a given s 

Ci{s) = -J^ dso{t{so) ■ t{so + s))^. (21) 

Hence, the zero-energy torsional fiuctuations associated with sliding the twist-wall ar- 
ray along the backbone distribute the weak spots uniformly along the chain and lead to 
translationally-invariant tangent correlations. 

In Fig. [3] (b) we plot angular excursion, (|6'(0, s)p)^, for filaments with a bending aspect 
ratio, C1/C2 = 16 and = 0, for weak and strong surface binding. Above the absorption 
threshold, we find that angular diffusion grows linearly with arc-length, with the relative 
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small slope of l/£i. Just below the critical binding strength, the signature of twist walls 
appear, where the local persistence length drops to £2 = ^i/4 over a relativity short length 
span of roughly A, corresponding to "weak spots" of easy bending around the thin direction. 
As the binding strength drops to negligible levels, the filament adopts its native twist, uq, and 
the local persistence length, as measured by the inverse slope of {\9{0, oscillates between 
ii and £2 over equivalent quarter-pitch spans of arc-length. Fig. |4] (a) shows — InC^(s) for 
a range of binding strengths, for the case iico = 1 and C1/C2 = 16 as in Fig. |3] (b). Despite 
the translation invariance of C^(s) regularly spaced intervals of relatively flexible backbone 
are still apparent in the tangent-tangent correlation functions below the threshold binding 
strength. 

Though the angular diffusion is not strictly a linear of function of arc-length over short 
distances, we may extract the long-distance, effective persistence length from the fall of 
tangent-tangent correlations over a span of P/2 corresponding to a period span of the arc- 
length modulation, 

{mp/2r)jp. (22) 



lnCt(P/2) _ ^,,,,2, 



In Fig. |5] we plot the dependence of 4ff on the binding strength, for three values of 
cross-section anisotropy. Clearly, in the locked-in, untwisted state, the persistence length 
is maximal, £efr = ^1, as in-plane bending occurs only around wide axis of the filament. 
In the opposite limit of vanishing binding strength, the uniform helical rotation of the 
filament, ip = uqs, samples bending around all cross-sectional directions evenly, and from 



eq. (20) we find that ics = (^1^2)^^^- Hence, going from strong binding to weak binding 
of a helical filament we observe a (^1/^2)^^^ = w/t-fold drop in the persistence length of 
the backbone. From Fig. |5] we see that most of this drop occurs over a narrow region near 
to the untwisting transition, V ^ Vc- However, it is significant to note from Fig. [s] that 
even far from this transition, at weak binding strengths we predict a weaker, and roughly 
linear, dependence of 4ff on V. This demonstrates that the apparent persistence length of 
a surface absorbed helical filament is always sensitive to anisotropic effects of surface binding. 



We know look at actual conformations of the polymer as we cross this transition. Monte 
Carlo simulations revealed an excellent agreement with analytic results as seen in Fig. |3] and 
Fig. [5} Moreover the presence of twist walls (Fig. |4]) can be seen in our simulation results 
as well. In Fig. |6] we can see two chain conformations which characterize qualitatively the 
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FIG. 5: Plots of the effective persistence lengths as functions of surface binding strength for 
filaments of different in-plane anisotropy. For a continuum model of elasticity, a given width and 
thickness ratio implies a ratio of bending moduli, C\ = {'w/t)'^C2- Monte Carlo simulation data 
with vertical error bars are plotted on top of analytical result and show good agreement with 
theoretical predictions. 

interplay of average twist and effective persistence length. Fig. [6] (a) is a snapshot of a 
collapsed chain configuration when V" = 0. The chain twists at its natural twist rate and 
the effective persistence length is leading to a collapsed chain state. In contrast Fig. |6] (b) 
shows a situation when the binding potential is above the critical potential. Here the chain 
is almost fiat as predicted and the effective persistence length is large leading to a swollen 
chain configuration. 

To quantify the collapse transition in a more experientially accessible way, we now focus 
on the squared end to end displacement, {R'^))/L^ normalized by the contour length of the 
polymer. For binding potentials V < VcV/e can see that (i?^))/L^ increases as the potential 
is increased from zero. This behavior should be expected from the behavior of the ies data 
(Fig. [5]). Once V > Vc the polymer untwists and only the larger of the two bending stiffness 
is accessible when {ip') = and hence the fluctuations of the polymer should be indicative 
of that length scale. This picture is consistent with the data in Fig. [7] in the V > Vc regime 
where is constant. An expression for the squared end to end displacement for a 

fixed contour length in 2d is [36] , 
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FIG. 6: Typical Monte Carlo chain configurations acquired during simulation runtime at two values 
of binding potential, (a) For V = the chain twist with its natural twist rate loq and is collapsed 
with a characteristic persistence length given by eq. ( |12| ). (b) For V > Vc we see that twist walls 
are being expelled and the chain is essentially flat with ii as the characteristic persistence length. 
Here C1/C2 = 2 and L=100 subunits. 



(R^) 44, 2L, , X 

^-j/ = ^{l-^{l-e^ "''0) (23) 

In Fig. [7] the dashed lines were generated using the extracted persistence lengths from 



the data in Fig. pi As we can see the {E?) j L"^ data is in very good agreement with eq. (23). 



III. CONCLUSION 



Structural and mechanical properties of biopolymers adsorbed to a surface are of great 
importance in biology. In vivo there are many examples of membrane bound polymers 
which are essential for cellular functions such as cell shape maintenance, motility, mechanical 
sensing and the cleaving and contraction of the cell membrane during division. It is therefore 
imperative to understand if and how biopolymers change their mechanical and structural 
properties when adsorbed/bound to a membrane surface. This question becomes particularly 
interesting when one considers that real biopolymers typically have helical structures and can 
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FIG. 7: Simulation data for the end to end displacement normalized by the square of the polymer 
contour length. Each data set if for a fixed contour length but varying aspect ratio. As the aspect 
ratio of the polymer is increase from w/t = 2 (red) to w/t = 8 (blue) the magnitude of the 
transition to constant persistence length [V > Vc) is less dramatic. Dashed lines are generated 



have anisotropic bending stiffnesses. The interplay between the intrinsic helicity, the binding 
affinity to the surface and the bending stiffness anisotropy lead to non-trivial conformational 
and mechanical changes upon filament binding, which in turn can be tuned or regulated by 
the filament's affinity for the surface. 

Here we have studied a continuum model for a ribbon polymer which has an anisotropic 
bending stiffness and a native helical twist around it's long axis. While the anisotropy 
in a real biopolymer can originate in the detailed chemical structure, here we choose to 
model the biopolymer as uniformly elastic ribbon with a width w and thickness t, thus the 
structurally the anisotropic bending stiffness arises when aspect ratio w/t ^ 1. The helical 
twist implies that the filament's affinity for the surface will be modulated on a period set 
by the native twist rate when the filament is in its naturally twisted state. In the context 
of the ribbon model we assume that the wider side with more surface area also has a higher 



using the ^eff data from figure [5] and eq. ([23]). 
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binding affinity for the surface and the twist therefore determines which side is in contact 
with the surface. Though this assumption is follows intuitively from the fact that face-on 
binding of the wider edge of the cross-section allows for greater filament-surface contact, 
the detailed orientation dependence of ffiament-surface interactions need not conform to 
this simplified picture. It is quite easy to see that, however, that any sufficiently strong 
orientation dependence of interactions (say, a strong preference edge-on binding) will give 
rise to the unwinding transition described here. 

We ffist considered an infinitely stiff limit {9' 0) where the larger surface area is in 
contact with the binding surface most of the time. For surface potential strengths above 



the critical value Vc (eq. 17), the net twist vanishes and the wider face of the ribbon 
maintains contact with the surface. Once below the critical value of the potential {V ^ Vc), 
a sharp continuous transition occurs to a configuration where twist walls proliferate and 
form a periodic array (in the absence of thermal fluctuations) which are are separated by a 
distance P/2. Due to the highly heterogeneous structure of the twist in this transition from 
the weakly-bound to strongly-bound regime, we found that spectrum of in-plane bending 
fluctuations is highly sensitive to surface interactions near to this transition. In the weak- 
binding limit, thermal bending of the backbone is enhanced in the "softer" regions where 
the filament has an edge-on orientation. Here, excursions of the in-plane angle 6 oscillate 
around a purely diffusive dependence in the contour length, a signature of the appearance 
of twist walls which function as floppy joints with the period of the oscillations being set 
by the average distance between twist walls . The effective persistence length has a highly 
nonlinear increase with increased potential for V < Vc- As the magnitude of the potential 
approaches the critical potential from below we find that the amplitude of the oscillations 
diminishes to zero as the period continuously approaches infinity. Mechanically this is due 
to the larger fiat side constituting more and more of the bound side and thus the effective 
persistence length approaches that of the stiffer bending direction, leading to a homogeneous 
bending stiffness corresponding to the bending about the wide ffiament direction. 

We propose that one could use the strong dependence of the apparent persistence length 
on surface interactions to address in important and previously unresolved biomechanical 
property of biological filaments experimentally. Though the bending stiffness of helical 
bioffiaments like F-actin is well characterized on long lengthscales, the are no measurements 
to our knowledge of the local anisotropy of bending stiffness. By modulating the strength of 
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filament-surface interactions and monitoring fluctuations of the backbone orientation, say by 
modulating surface charge or osmotic pressure of the solution, the dependence of effective 
persistence length on binding strength, iefi{V) can be measured by standard microscopy 
techniques. Extrapolating these data to the limits of low- and high-surface binding one has 
a direct access to stiffness anisotropy via the relation, 4fr(^ -> oo)/4ff(^ ^ 0) = ^Ci/C2. 

To assess the physiological relevance and experimental accessibility of surface-induced 
untwisting helical biopolymers, we now discuss known values of the physical parameters 
affecting this behavior. For simplicity we focus on electrostatic interactions between the 
polymer and surface, which is easily accessible in vitro. We are interested then in the 
surface charge density, a, required to attain the critical binding strength Vc say for F-actin 
whose long-wavelength mechanical properties are well-characterized. F-actin is negatively 
charged with linear charge density of p ~ 4 nm~^ (in units of e) [37], and thus in solution 
would bind to an oppositely charged surface. Given the equation for the critical binding 



strength (eqn, 17) we can calculate the necessary charge density. First we quantify the critical 
potential per unit-length that is required to untwist F-actin, we use a native twist of 1.6 
turns per micron (10 fim~^) [3H] and a twist stiffness in the range of 8 x 10^^ pN /zm^ This 
provides us with a scale of critical potential (per-unit length). 



K ~ 10 pN (24) 

To quantify the critical surface charge we assume that the solution of F-actin is dilute 
such that the distance between neighboring adsorbed filaments is very large compared to 
their length. Assuming a screened electrostatics, the electrostatic free energy gained per 
unit charge brought to an oppositely charged surface is roughly fesTcr^s//?, where is is the 
Bjerrum length and k,~^ is the Debye screening length of the solution. From this we may 
crudely estimate the value of the surface-binding potential to be of order 

V keT ^ . 

From this we estimate the value of cr for which we expect surface interactions to be 
sufficiently strong to drive the unwinding transtion of F-actin, a ~ Vcn/ksTpiB, at room 
temperature in 1 fiM monovalent salt, where is/i^ = 0.21 nm^. The critical surface charge 
to fully untwist the filament is, 
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ffc = 3.e nm 

These values are consistent with physiological values for the surface charge density on 
cell membranes indicating that such conformational transitions and the associated changes 
in mechanical properties of surface bound filaments may in fact be exploited in vivo. The 
modest scale cTc required also suggest that such the a careful regulation of electrostatically- 
induced surface binding can be used to indirectly probe the anisotropic mechanical properties 
of helical helical biofilaments like F-actin and bacterial homologs like MreB. It is important 
to take these effects into account when backing out elastic properties of twisted filaments 
adsorbed on surfaces. For instance, there have been recent observations that DNA can 
form abnormally high bending deformations when adsorbed to a surface [39]. DNA is also 
intrinsically twisted (1 revolution per 3.4 nm) on length scales where anomalous elasticity 
is observed (~ 5 nm). In the context of our model we have presented here it would be 
interesting to understand how local denatured "weak spots" could affect the formation of 
twist walls and how the overall observed persistence length is determined by the interplay 
of intrinsic twist and adsorption onto a surface. 
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Appendix A: Twist profile solutions 

Given the constant of integration, e, the relationship between ijj and arc-length follows 



the solution to eq. (13) 



s = X r d^' I = Xe-'/'Fi^P, -e-i), (Al) 

Jo V sin yj' + e 
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where k) is the incomplete elhptic integral of the first kind and we have set ipls = 0) = 0. 
This equation can be inverted to given ip{s) in terms of the Jacobi elliptic amplitude function, 

's^/e 



ipi-s) = am 



A 



— e 



-1 



(A2) 



Using the solution to Euler-Lagrange equation for ipi^), "we may derive an expression for the 
energy of a single pitch of the bound helix as function of the parameter e, 



E( 



P{e) 



ds 



K. 



(A3) 



Using eq. (13) we may the integral above as 



ds 



IV- I +1^sin^^ 



2-K 

2V\ I dp\/sm'p + e-VeP{e) 





From this, we have the energy density, 

E{e) _ 2V\ '■^^ 



dip A/ sin ip + t — Vt 







Pie] 



(A4) 



(A5) 



which we minimize with respect to e to derive the equation of state, effectively relating mean 
pitch to the binding potential and twist-elastic cost for V <Vc, 



ttKuq 
VX 



[ dipJsin^ij + e = 4e^/^E(-e~^), 
Jo 



(A6) 



where E{k) is the complete elliptic function of the second kind. For a given e this relation 
can easily be solved to binding strength, 

Vr. 



Vie) 



(A7) 



Noting that A = XdV/Vcy^/^, where A^^ = tiujo/2, eqs. (A7) we may rewrite the equilibrium 
pitch in terms of complete elliptic integrals. 



Pie) = KEi-e'')Ki-e 



(AJ 



Using these expressions. Vie) and -P(e), we describe the relationship between equilibrium 
pitch of the helical filament and binding strength through the parametric dependence on e, 
plotted in Fig. [31(b). 
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